Reply to the comment on: 
"Thermostatistics of Overdamped Motion of Interacting Particles" 

by Y. Levin and R. Pakter 
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In their comment, Levin and Pakter present an ana- 
lytical solution for the mechanical equilibrium of a sys- 
tem of particles interacting through a potential given by 
the modified Bessel function, and confined by a restor- 
ing force. This derivation yields a quadratic dependence 
for the local density of particles, which is consistent with 
Eq. (14) of our work. Based on this result, Levin and 
Pakter question our interpretation that a system of over- 
damped particles at T=0 is a physical realization of Tsal- 
lis thermostatics with entropic index v = 2. They claim 
that the same density profile can be obtained using only 
Newton's laws. 

In what follows we provide a reply to the remarks of 
Levin and Pakter. We first show that their attack to 
our results and analysis is conceptually unfounded and 
rather misleading. Inexplicably, they simply choose to 
categorically dismiss our elaborated and solid concep- 
tual approach and results, without employing any fun- 
damental concepts or tools from Statistical Physics. We 
then demonstrate that the results of Levin and Pakter 
do not present any evidence against, but rather corrob- 
orates, our conclusions. In fact, the results shown in 
their comment correspond to a confining potential that 
is 1000 times stronger than the typical valued utilized in 
our study, therefore explaining the discrepancy between 
their results and ours. Furthermore, in this regime where 
higher vortex densities are involved, vortex cores might 
get so close to each other that can no longer be treated as 
point-like defects. As a consequence, Ginzburg-Landau 
equations should be employed instead, meaning that the 
physical conditions implied by the results of Levin and 
Pakter should be considered with caution in the context 
of the Physics of interacting superconducting vortexes. 

1) Levin and Pakter claim that our result for 
the density distribution of particles "has nothing 
to do with the Tsallis statistics" . 

In our letter, we show that the the non-linear Fokker- 
Planck equation in the form 
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results in a dynamics that drives the system towards a 
state that maximizes Tsallis entropy with entropic index 



v. See also [l[ where this is discussed in more detail. Note 
that Levin and Pakter do not contest this statement. 

We also show that, for a class of systems of overdamped 
particles interacting through repulsive short range poten- 
tials, the particle-particle interaction induces a flux of 
particles proportional to gradient of concentration. This 
has been previously discussed in few other works Q • We 
also show in our letter that this kind of dynamics leads 
to a non-linear diffusion equation that is equivalent to 
the non-linear Fokker-Planck equation ([lj. Again, none 
of these statements were contested by Levin and Pakter. 

In view of the aforementioned facts we are rather sur- 
prised that Levin and Pakter, without contesting any of 
the results demonstrated in our work, still question our 
conclusion that "a system of overdamped particles at T=0 
is a physical realization of Tsallis thermostatics with en- 
tropic index v — 2." This is specially puzzling when 
one notes that the solution proposed by Levin and Pak- 
ter for this system also corresponds to a solution of our 
non-linear equation, therefore confirming our findings. 

2) Levin and Pakter state that our result for 
the density profile of particles "has everything to do 
with Newton's Second Laws" . 

Of course, we were by no means expecting our results 
to be in contradiction with Newton's Laws. As a mat- 
ter of fact, in our study, we used molecular dynamics 
to integrate equations of motion and obtain numerically 
density profiles, validating the predictions of our non- 
linear Fokker-Planck equation. In the same way, neither 
the Maxwell-Boltzmann nor the Tsallis thermostatistics 
are contrary to classical mechanics theory, as one could 
infer from the final and rather misleading sentence of the 
comment. 

3) In their comment, Levin and Pakter also 
state that "prior to discarding the standard statistical 
mechanics one should see what it has to say on this mat- 
ter:' 

Again, our approach and results neither disregard nor 
contradict "standard" statistical mechanics. We would 
surely welcome any other treatment based on the solid 
principles of Statistical Physics. This is definitely not the 
case of this comment. Instead, however, they categori- 
cally dismiss our elaborated approach and results without 
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employing any concept or tool of Statistical Physics. 

4) Levin and Pakter write in their comment 
that "if a classical system is placed in contact with 
a temperature reservoir at T = it will loose all its 
kinetic energy and collapse to the ground state" 
where the "net force on each particle vanishes." 
They also state in their comment that "this is 
precisely what happens for the system studied by 
Andrade et al." 

Of course, the overdamped motion we investigate 
should evolve towards mechanical equilibrium, as we ex- 
plicitly state in our paper. One should note, however, 
that our approach goes beyond the prediction of a sta- 
tionary state, since the non-linear Fokker-Planck equa- 
tion introduced by us also describes the dynamics of the 
system, with all its transient features, till the station- 
ary regime is eventually reached. The analytical solution 
provided in the comment, however, does not contemplate 
the dynamics outside the stationary state. The fact that 
our system evolves towards mechanical equilibrium is in 
no way contradictory to our results or conclusions. 

It is important to point out that we also studied the 
case in which thermal noise disrupts equilibrium, result- 
ing in a density function that is consistent with a mixed 
entropy combining Tsallis with entropic index v = 2, and 
Boltzmann-Gibbs formulations, namely, Eq. (17) of our 
paper. By varying the ratio between the particle-particle 
interaction and thermal force, we can continuously inter- 
polate the density profile from a Gaussian, at higher tem- 
peratures, to a parabolic shape at T — 0. In short, what 
we proposed in our paper and the authors of the com- 
ment appeared to miss completely, is that overdamped 
particles interacting through a wide family of potentials 
follow a dynamics that drives the system to the maximum 
value of the entropic function Eq. (17) of our paper. The 
usefulness of our approach is that the density profile can 
be easily obtained with the many-body interactions being 
carried into the entropic term with v = 2. 

5) In their comment Levin and Pakter claim 
that they obtain an exact solution for the particle 
distribution at T = 0. 

This claim is misleading to say the least. In fact, in 
the same way as we did, their solution is a coarse-grain 
approximation, where the discrete nature of the interact- 
ing particles is replaced by a continuum description for 
the particle density. However, it is relevant to note that 
they arrive at a solution that is similar, but not identical 
to ours. Levin and Pakter also show results from numer- 
ical simulations that seem to agree with their solution. 
As we show here, however, they studied the system at 
a regime that is significantly different from the one we 
investigated in our work, and that is clearly the reason 
for the discrepancies. 

The equation of motion for our system of overdamped 
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FIG. 1. Numerical results for the density profile obtained 
from numerical simulations. We used the same conditions as 
Fig. 1 of our letter, namely, N = 800,q = 10~ 3 / , L y = 20A. 
The dashed red line is the solution for this profile given by 
Eq. (14) of our letter, while the straight blue line correspond 
to the solution according to Eq. (3) from Levin and Pakter's 
comment. Clearly their solution is not compatible with the 
concavity of the density profile obtained from our molecular 
dynamics simulations. The inset shows that the discontinuity 
in their solution is barely noticeable for these physical condi- 
tions. 

particles is given by 

M = / ] foKiQfjj 1/ X)fjj - axj (2) 

jjti i 

where K\ is the modified Bessel function. The factors 
/o and a control the intensity of the particle-particle in- 
teraction and the confining force, respectively. In their 
results Levin and Pakter used a = fo while we used 
a = 10~ 3 /o. It is not surprising that, by using a con- 
fining potential 1000 times stronger, one should obtain 
results that are distinct from ours. 

6) Levin and Pakter arrive at a result that is 
similar to Eq. (14) of our paper, namely, Eq. (3) of 
their comment. They claim, however, that their 
solution does not need any adjustable parameter, 
while we had to set the value of the parameter 
that controls the concavity of the distribution to 
a = 2.41/ A, in order to adequately fit our model to 
our results from molecular dynamics simulations. 

Apparently the authors of the comment did not realize 
that their proposal of a = tt/qX 3 as a "non-adjustable pa- 
rameter" corresponds exactly to the same theoretical pre- 
diction presented in our paper. However, this prediction 
is valid only when the density of particles p varies slowly 
within the interaction range of the potential. Since the 
system is composed of discrete particles, the value of the 
parameter a should approach 7r/oA 3 as the density grows. 
Since Levin and Pakter use a substantially stronger con- 
fining potential, the obtained densities can reach 5 par- 
ticles per A 2 , which is about ten times higher than ours. 
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At this regime, the value a = 7r/ A 3 represents indeed 
a good fit to the numerical data. At moderates densi- 
ties, however, the numerical results do not follow this 
prediction. In Fig. [TJ we show our numerical results and 
parabolic predictions with a = 2.41/ A 3 and a = 7r/ A 3 . 
It is clear that the value of a used by Levin and Pak- 
ter do not agree with numerical results from molecular 
dynamics in this regime. 

In any case, it is important to emphasize that the inter- 
acting potential in the form of the modified Bessel func- 
tion is originally motivated by applications in the theory 
of superconducting vortexes. It is a known fact Q for 
this physical system that such a potential form repre- 
sents an approximated model that is only valid in the 
regime of moderate magnetic fields (i.e., moderate par- 
ticle densities, p <C k 2 /2it\ 2 , where k is the Ginzburg- 
Landau parameter). In the case of higher densities, the 
vortex cores might get so close to each other that can 
no longer be treated as point-like defects. As a conse- 
quence, Ginzburg-Landau equations should be employed 
instead, meaning that the physical conditions implied by 
the results of Levin and Pakter should be considered with 
caution in the context of the Physics of interacting su- 
perconducting vortexes. 

7) Levin and Pakter predict the existence of 
discontinuities at the edges of the density function 
P- 

We would like to point out that, although the numer- 



ical results of Levin and Pakter follow their theoretical 
prediction in the bulk of the density profile, by just look- 
ing at their histogram it is not clear whether the dis- 
continuity is really present of not. The inset of Fig. Q] 
shows that under the conditions of our simulations, this 
discontinuity, if present, would be negligible. Therefore, 
we conclude that, in the conditions we model our system, 
the discontinuity predicted by Levin and Pakter does not 
represent a relevant effect. In any case, the analytical re- 
sult of Levin and Pakter can be viewed as a solution of 
Eq. (13) of our letter, with a — 7r/ A 3 , but subjected 
to a different boundary condition. Therefore, our under- 
standing is that their results are also consistent with our 
theoretical model. 

In summary, based on the preceding comments and 
remarks, it is our opinion that the comment by Y. Levin 
and R. Pakter is neither conceptually sound nor relevant. 
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